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The properties of neutral and anionic Pdjv clusters were investigated with spin-density-functional 
calculations. The ground state structures are three-dimensional for 7V>3 and they are magnetic with 
a spin-triplet for 2<N<7 and a spin nonet for 7V=13 neutral clusters. Structural- and spin-isomers 
were determined and an anomalous increase of the magnetic moment with temperature is predicted 
for a Pd7 ensemble. Vertical electron detachment and ionization energies were calculated and the 
former agree well with measured values for Pd^. 



Enhancement of magnetism in finite clusters made 
of elements that are ferromagnetic as bulk solids has 
been demonstrated through Stern-Gerlach (SG) deflec- 
tion measurements Q, and it is understood to derive 
from reduced atomic coordination resulting in stronger 
electron localization. However, the emergence of mag- 
netism in small clusters of close-shell non-magnetic atoms 
is clouded by uncertainty. Palladium aggregates (with a 
[Kr] 4d 10 structure of the atom) are particularly interest- 
ing since the bulk metal is known to be "almost" mag- 
netic, with the emergence of magnetism predicted to re- 
quire a mere 6 % dilation of the interatomic distance . 
Experimental information is limited to SG measurements 
on Pdjv clusters with iV>12 (where the clusters were 
found to be non-magnetic) || or to inferences from (non- 
magnetic) photoemission studies Q|. Additionally, sys- 
tematic theoretical investigations of magnetic properties 
are lacking, limited only to neutral Pdi 3 ||||. 

We report on an extensive density-functional-theory 
(DFT) study pertaining to size-dependent evolutionary 
patterns of the properties of Pd^ and Pd^y clusters with 
N=l-7 and 2V=13. The ground states (GSs) of the neu- 
tral and anionic clusters are found to undergo an early 
transition (i.e. for 7V>3) to three-dimensional (3D) ones 
and possess a non-zero magnetic moment (e.g., for neu- 
tral clusters with N<7 a spin triplet, 5=1, and a spin 
nonet, 5=4, for iV=13). In addition to higher-energy 
structural isomers (STIs) we determined sequences of 
close- lying spin isomers (SPIs). For a Pd7 ensemble, we 
predict an increase of the magnetic moment with tem- 
perature due to the thermal accessibility of such iso- 
mers. All the SPIs exhibit high local magnetic moments 
(LMMs), including the singlet states with an antifer- 
romagnetic LMM coupling. Our structural determina- 
tions are corroborated by the remarkable agreement be- 
tween the calculated vertical electron detachment ener- 
gies (vDEs) from the cluster anions and photoelectron 
spectroscopy (PES) measurements 

In this study the Kohn-Sham (KS) equations with 
generalized gradient corrections (GGA) M were solved 



using the Born-Oppenheimer local-spin-dcnsity molecu- 
lar dynamics (BO-LSD-MD) method pQ |, with scalar- 
relativistic |ll]] non-local pseudopotentials |Q. Cluster 
geometries were determined via symmetry-unrestricted 
structural optimizations through a rather exhaustive 
search starting from various structures including those 
suggested in previous studies of certain Pd clusters 
jl3] jl5| and other metal clusters 16 17|]. For each struc- 
ture spin-restricted optimizations were performed cover- 
ing all energetically important spin multiplicities. 

The GS geometries and STIs for both Pd^v and Pd^y 
are rather similar and they follow the same structural 
evolution, exhibiting Jahn- Teller distortions from the 
ideal symmetric structures (see Fig. 1), with a transi- 
tion to 3D configurations at iV>3 [§J|]. The GSs (i.e. 
maximal binding energy Eb in Fig. 1) of Pd^ with 
N<7 have a triplet spin multiplicity (see open squares 
in Fig. l(A-F)), while for Pdi3 the GS is associated with 
a nonet (5=4) spin configuration (see Fig. 1G). On the 
other hand, the GS spin-multiplicities of the Pd^ an- 
ions (filled symbols in Fig. 1) vary non-monotonically 
with N (doublet for 7V=2,3,4, an d 6; quartet for A=5; 
sextet for N=7, and octet for iV=13). We note here 
that the higher-lying STIs |lj| and SPIs of the neutral 
and anionic clusters become thermally accessible with in- 
creasing cluster size (see temperature scales |2TJ on the 
right-hand-side in Fig. 1). 

The binding energies of the GS clusters increase rather 
monotonically with N (Fig. 2a), showing enhanced local 
stabilities for Pdg and Pd 4 The average nearest- 

neighbor bond- length ((R nn )) of the GS structures and 
the SPIs are very similar [p2[ , converging rapidly to the 
bulk value (2.75 A @). 

The calculated vDEs from the cluster anions are in a 
remarkable agreement with values determined from PES 
measurements (Fig. 2b). For all the clusters, we 

display the vDEs only for the GS (structural and spin) 
configurations of the anions, except for Pdy where we 
show the vDEs for both the GS (5 = 5/2) and its next 
(higher- in-energy) SPI (5 = 3/2) (the difference in en- 
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ergy between the two isomers is 0.11 eV, corresponding to 
only 170 K in vibrational temperature p(J), with the lat- 
ter exhibiting a better agreement with the experimental 
data. These results suggest that for Pdf the measured 
vDE may correspond to this higher-energy SPI which is 
accessible already at very low temperatures, while for the 
other clusters studied here the vDEs are determined by 
the GS structures over a broad temperature range. 

The only (slight) discrepancy between the calculated 
and measured VDEs is for Pd^ (see Fig. 2b). How- 
ever, we observe from this figure that the predicted vDE 
(2.0 eV) for Pdg~ is rather close to the dissociation en- 
ergy (Ed) of the process Pd^— >Pd^~+Pd (see filled circle 
in Fig. 2b at Ed=2A\ eV, compared to a measured 
value of 2.26 eV 63). This implies that this dissociation 
channel may compete with the electron detachment un- 
der appropriate experimental conditions. Consequently, 
we assign the vDE of Pdg to the next higher feature in 
the measured PES ||) occuring at 1.88 eV which agrees 
with the predicted value, and the lower measured value 
(at 1.66 eV, marked by a cross in Fig. 2b) is attributed 
to the vDE of a (hot) Pd^~ dissociation product (the 
same conclusions apply to the photoelectron spectrum of 
Pd^ shown in Fig. 2 of ref. Q); for a similar interpre- 
tation of the measured PES of Aug" see Refs. 16 and 
26. 

The vertical ionization potentials predicted by us (vIP 
in Fig. 2c) start at the atomic value of 8.28 eV com- 
pared to the measured first IP of Pd, 8.33 eV |23|; the 
calculated second IP of Pd (27.34 eV) also agrees with 
the experimental one (27.75 eV) p3| . The predicted val- 
ues for N—2 and 3 are close to each other followed by 
a marked drop for clusters with 7V>3, which is likely to 
reflect the transition to 3D structures. Convergence to 
the bulk limit (the work function of Pd is 4.97 eV Hf)) 
is slow. No measured IP values have been reported. 

The multitude of spin-multiplicities for the Pd^ clus- 
ters results in a considerable variation of the magnetic 
moment per atom /i=25 /is I 'N (solid squares in Fig. 2d) 
with particularly high values of /j,=0.6, 0.71 and 0.54 
for N=<3, 7 and 13, respectively. On the other hand, the 
occurrence of a triplet GS for the neutral clusters under- 
lies a monotonic 1/N decrease of \i for N=l — 7 (open 
squares in Fig. 2d), with an unexpected high value of 
0.62 fig for Pdi3 which is higher than the experimen- 
tal H estimate (/i<0.4 fis)- 

The rather surprising existence of energetically favor- 
able high-spin isomers for Pd clusters originates from 
sd-hybridization developing upon bonding, with the d- 
orbitals' weight diminished somewhat on each atom, re- 
sulting in a situation reminiscent of open-shell transition 
metals. The total magnetic moment of the cluster is com- 
prised of sizable atomic LMMs (/i^~±(0.3 — 0.6)hb) that 
couple antiferromagnetically in the spin-compensated 
singlet states, and align themselves in SPIs with high 
total fi. This is illustrated for the icosahedral Pdi3 clus- 



ter in Fig. 2 where the (gapless) density of states (DOS) 
of the 5 = singlet SPI is shown in Fig. 2e and the 
corresponding spin-polarization density in Fig. 2g (left); 
note the non-uniform spatial distribution of the spin- 
polarization reflected in the different line-shapes of the 
DOS of the up and down spins. In the 5 = 2 spin-quintet 
SPI the minority spin polarization is localized on three 
sites located in a triangle (Fig. 2g, middle). In the 5 = 4 
spin-nonet GS cluster all the sites are spin-polarized in 
the same direction (Fig. 2g, right), and it's stability is 
reflected in the large gap in the majority-spin DOS near 
the Fermi energy (see spin f in Fig. 2f). 

Intriguing conclusions can be made regarding the ther- 
mal behavior of certain Pd clusters in SG measurements; 
here we consider the case of Pd7 and Pdy at room tem- 
perature. Neglecting the vibrational and entropic differ- 
ences between the isomers, the probability to find a clus- 
ter with spin 5, irrespective of its atomic isomeric struc- 
ture, is given by P n ,t(S) = Y,i ex -V(^^r J ^-)l z N.T 
with the Boltzmann constant ks, the ensemble temper- 
ature T, the structural isomer index J, and the (normal- 
izing) partition function Zn,t- Several spin isomers of 
the neutral as well as the anionic heptamer have a finite 
P?,t(S) (Fig. 3a) even for low temperatures, leading us 
to predict that in a SG experiment up to 3 different de- 
flection angles should be measurable. Note, that while an 
increase in temperature depopulates the Pdy pentagonal 
bipyramid triplet (5=1) GS, the singlet (5=0) state of 
that structure and the capped octahedron quintet (5=2) 
state (Fig. IF) gain statistical weight. Consequently, we 
predict that for Pd7 a rise in temperature would lead 
first to a decrease of the thermally averaged magnetic 
moment per atom (h)n,t=J2s HbPn,t(S) /N due to a 
sharp increase in the population of the singlet state. The 
subsequent increase of (h)n,t (see Fig. 3b, T>200 K) 
results from the higher thermal population of the quin- 
tet state relative to the GS triplet (Fig. 3a). Q Such 
anomaly does not occur in the case of the Pd? cluster 
where the doublet (5=1/2) and quartet (5=3/2) states 
start to coexist with the sextet (5=5/2) GS at elevated 
temperatures (Fig. 3a). For both neutral and anionic 
Pdi3, the higher- lying isomers play essentially no role for 
T<800 K resulting in a weak temperature dependence of 
(h)n,t (Fig. 3b). 

In summary, we found that unlike atomic and bulk 
Pd, both neutral and anionic Pdjv clusters (2 < N < 7 
and N = 13) are magnetic, with relatively high LMMs. 
Underlying this behavior is the hybridization of atomic 
s and d states when clusters are formed, that depletes 
local d-contribution around each atom and leads to an 
open-shell-like behavior. The abundance of close-lying 
SPIs for certain clusters should be detectable in ther- 
mally controlled SG experiments, and we predict an en- 
semble of Pdy clusters to exhibit an (anomalous) in- 
crease of the magnetic moment with temperature. The 
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remarkable agreement between the calculated and mea- 
sured vDEs of the cluster anions corroborates the pre- 
dicted atomic structures. Our results provide the first 
quantitative predictions pertaining to the emergence of 
pronounced magnetic properties of Pdjv and Pd^ clus- 
ters, unlike the weak magnetic tendencies inferred indi- 
rectly from PES data 0] ; note in particular the contrast 
between our results for Pd7 (see Fig. 2d and Fig. 3) 
and the suggestion of a zero spin GS for Pd7 given in 
Ref. 01. Furthermore, our study motivates temperature- 
dependent magnetic deflection (SG) measurements and 
further investigations of free and supported Pd clusters, 
including correlations between their magnetic properties 
and their catalytic activity |^0| . 
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FIG. 1. Binding energies Eb (in eV) of structural (a,b,c) and spin (5=0, 1/2, 1, 3/2 etc.) isomers of Pd^ and Pdjv for 
N=l~7 (A-F) and 13 (G). Filled and open symbols represent anions and neutrals, respectively with squares for the ground state 
structure (a), circles for structure (b) and triangles for (c). The anionic Eb is defined as the energy per atom to separate Pd^ 
into -/V neutral atoms and an electron. The Eb values on the left hand side of the panels are mapped on separate temperature 
scales |2(j] (T, in K) for the anions (top values on the right hand axis of each panel) and the neutrals (bottom values), giving 
an estimate of the thermal accessibility of the isomers. Note the large number of isomers for Pd7 and Pdy for T<400 K. Some 
of the structural isomers were not stable for certain spin values (such cases are denoted by dashed lines). The PdJ square (b) 
in panel C transformed into a rhombus (c) for 5=3/2 and the Pd4 rhombus transformed into a square for 5=0 and 1. The 
Pd6 capped trigonal bipyramid (b) transformed to the octahedron (a) for 5=0. The Pd7 capped octahedron (b) became a 
pentagonal bipyramid (a) for 5=3 and the Pdi3 and Pdj~ 3 cuboctahedron (b) transformed to the icosahedron (a) for 5=4.5 
and 5, respectively. 

FIG. 2. Size evolution (panels a-d) of energetic and magnetic quantities in Pd clusters, (a): binding energy per atom, Eb 
(in eV), (b): theoretical (filled squares) and experimental (open circles) vertical electron detachment energy (vDE) and the 
calculated atom dissociation energy Ed of the anions (filled circles). Experimental values are take from ref. |7U8| for 7V=1— 3 
(the corresponding results shown for PdJ in Fig. 2 of ref. [^] are essentially the same) and from Q for N—4-7 and 13. The cross 
and the open circle at JV=3 correspond to the maximum of peak A and the peak of group B-H in Fig. 1 of ref. 0, respectively 
(and correspondingly for the first and second peaks in Fig. 2 of ref. Q). For Pd^T, we show the vDE of the GS (5=5/2) and a 
thermally accessible SPI (5=3/2). The vDE of Pd^ was estimated according to Janak's theorem [^8| by the HOMO energy of 
Pd 1 1 , (c): calculated vertical ionization potentials; (d): the magnetic moment per atom of the GS anions (filled squares) and 
of the GS neutrals (open squares). (e,f):Density of states (DOS, in 1/eV) of singlet (e) and nonet (f) icosahedral Pdi3 cluster, 
(g): Constant-value images of spin-polarization density in singlet, quintet, and nonet icosahedral Pdi3 (left, middle, right, 
respectively) cluster. The purple and yellow denote excess of minority and majority spin, respectively. Note the transition from 
antiferromagnetic to ferromagnetic ordering when going from the singlet to the nonet. 

FIG. 3. (a): population probabilities of the different SPIs (irrespective the structure) of neutral (dashed curves) and anionic 
(solid curves) heptamers; (b): thermally averaged magnetic moment per atom for Pd7, Pdf, Pdi3 and Pd^ 3 . 
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